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Abstract 



The calculation of the hindered roton-phonon energy levels of a hydrogen 



' molecule in a confining potential with different symmetries is systematized for 

"j^ ■ the case when the rotational angular momentum J is a good quantum number. 

One goal of this program is to interpret the energy-resolved neutron time of 
' flight spectrum previously obtained for H2C6o- This spectrum gives direct infor- 

o ■ 

O ■ mation on the energy level spectrum of H2 molecules confined to the octahedral 

. interstitial sites of solid Ceo- We treat this problem of coupled translational 

^ , and orientational degrees of freedom a) by construction of an effective Hamilto- 

nian to describe the splitting of the manifold of states characterized by a given 
value of J and having a fixed total number of phonon excitations, b) by nu- 
merical solutions of the coupled translation-rotation problem on a discrete mesh 
of points in position space, and c) by a group theoretical symmetry analysis. 
Results obtained from these three different approaches are mutually consistent. 
The results of our calculations explain several hitherto uninterpreted aspects of 
the experimental observations, but show that a truly satisfactory orientational 
potential for the interaction of an H2 molecule with a surrounding array of C 
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atoms has not yet been developed. 



lypeset using REVI^ 



I. INTRODUCTION 



The study of rotational and vibrational dynamics of guest molecules (i.e. CO, O2, H2, 
etc) trapped in porous media such as fuUerenes, zeolites, and graphite has recently become an 
active subject both experimentally and theoretically.^"^ This is because such studies can yield 
valuable information about the host-guest interactions which could be important for several 
technical applications such as gas separation and hydrogen storage. In particular hydrogen 
molecules trapped in interstitial cavities in solid Ceo as well as hydrogen molecules embedded 
in nanotube ropes are of interest due to quantum behavior of hydrogen molecules in quasi zero 
and one dimensional sites. ^'^'^ 

In this paper, we develop a detailed analysis of coupled rotational and vibrational dynam- 
ics of a molecular hydrogen encapsulated in a solid using numerical, perturbative, and group 
theoretical methods. In particular we will be interested in what one might call the "weak cou- 
pling limit," when the interaction between molecular rotations and center-of-mass translations 
is weak enough that the rotational angular momentum quantum number J is a good quantum 
number. This limit is almost never satisfied except for very light molecules like hydrogen or 
deuterium. The energy levels of a free rotator are 

Ej^BJ{J + l), (1) 

where B — h^/(2I), where / is the moment of inertia of the molecule and Ej is (2 J -|- l)-fold 
degenerate. For H2 the rotational constant 2B has the vahic 60 cm~^, 14.7 mcV, or B/k = 85 
K (and the corresponding values for D2 are half as large), so that the energy separation between 
different J levels is large enough that often J is a good quantum number. This is certainly true 
for sohds consisting of these molecules unless the pressure is quite large. (For a review of the 
properties of the hydrogen molecule and sohd hydrogen see Ref. 6.) 

We have been led to consider this phenomenon in view of an experimental study of en- 
ergy spectra of H2 and D2 inserted into the octahedral interstitial sites in solid Ceo carried 
out by neutron time of flight techniques.^ In considering this phenomenon we should keep in 
mind the following experimental facts concerning the host solid of Ceo- The centers of the 
Ceo molecules form an fee lattice.^ At temperature above Tc, where Tc is about 260 K, the 
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molecules are orientationally disordered. At T = T(, long range oricntational ordering occurs^ 
and the molecules are ordered into four sublattices as described by Pa3 symmetry. In the 
orientationally disordered phase the local symmetry at the the octahedral interstitial site is 
indeed that of the point group Oh- In the presence of orientational ordering the symmetry of 
what was the "octahedral" interstitial site is now reduced to a uniaxial symmetry, specifically 
that of point group Sq}^ In experiments, hydrogen molecules are stable in the octahedral in- 
terstitial site only for temperatures well below Tc (where the interstitial site does not actually 
have octahedral symmetry). 

While a general understanding of the time-of-fiight experiments was presented,^ some of 
the finer details of the experiment remained unexplained. For instance, the shift in the energy 
associated with ortho-para {J = 1 ^ J — 0) conversion in the interstitial relative to its 
value for free molecules was not understood. Also the feature in the energy gain spectrum at 
about twice the ortho-para conversion energy was not unambiguously identified. These issues 
are both addressed in this paper. More generally we give a calculation of the cross section 
for neutron energy loss for comparison with the observed time-of-flight spectrum. For that 
purpose we need not only to consider the cross section for para-ortho conversion as compared 
to phonon creation, but also to calculate the phonon excitations of {J — 1) molecules. These 
calculations require us to develop and implement a scheme for treating coupled translations and 
rotations. In this paper we present a systematic analysis of the simplest case of this coupling 
which occurs when the quantum number J characterizing free rotation remains a good quantum 
number. In that case the well-known numerical schemes for solving the translation problem 
can be easily extended to include the effect of the coupling to rotations. In addition, we also 
give analytic expressions obtained by treating this coupling within perturbation theory. As we 
will see, this analj^ic development enables us to interpret many of the numerical results in a 
meaningful way. In addition, we analyze in detail various simplified models which illustrate our 
group theoretical analysis of the symmetry present in the system of coupled translations and 
rotations. This analysis indicates that arguments for the degeneracies of coupled translation- 
rotation modes based on simple classical concepts are incorrect. In summary: in this paper we 
present an analysis based on numerical, perturbative, and group theoretical methods. 
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Up to now there have not been many theoretical studies of energy levels in such irregular 
geometries like the octahedral interstitial sites in Cgo- A notable exception is the work of van 
der Avoird and collaborators^^ on CO in Ceo- That work examined an even more complicated 
situation in which the rotational and translation degrees of freedom interacted strongly. As a 
result, the problem was analyzed numerically. In contrast, for the present problem FitzGerald 
et al.^ applied a number of analytic and semi-analytic techniques to the theoretical study 
the spectra of hydrogen molecules in Cgo- This paper may be regarded as an extension and 
systematization of their approach. 

II. GENERAL FORMULATION 

Clearly the first step is to establish a satisfactory potential for the intercalated hydrogen 
molecule. This potential function V{r, fl) gives the energy of a hydrogen molecule whose center 
of mass is at r and whose orientation is specified by Q = {9, 0). A convenient starting point is 
to use an atom-atom potentiaP^ to describe the interaction between each of the two hydrogen 
atoms and the atoms in the confining structure. Unless otherwise indicated, all the results 
reported in this paper are obtained from the same WS77 potential, —A/r^ + 5exp(— Cr), 
that is used in Ref. 5 (where A = 5.94 eVA^ B = 678.2 eV, and C = 3.671"^). 

In this paper we will mainly consider the octahedral interstitial site in solid Ceo, but many of 
the considerations apply with slight modification to molecules confined within other structures 
such as single wall carbon nanotubes.-*^^ The determination of the potential V{r, il) for H2 in 
solid Cao is discussed in Appendix A. From the numerical evaluation of this potential we have 
extracted the expansion coefficients when it is written in the following canonical form: 

F(r,0) = H(r)+ E E . (2) 

Z=2,4,... m=-l 

We assume (and it is generally true) that the orientational energies which are relevant are much 
less than the smallest energy difference between successive J levels of a molecule {10 B for an 
ortho molecule and 6B for a para molecule). Accordingly, we may consider only that part of 
the potential which is diagonal in J. When the potential is written in the form of Eq. (2), it 
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is easy to implement the truncation to terms diagonal in J. So for a fixed value of J we have 
the Hamiltonian Tij as 

\J^,){J^,\Yr''mJ^^'){JA 

= -^V^ + H(r) + BJ{J + 1) + Vjiv) , (3) 

where Vj{r) is the orientationally dependent part of the potential (the terms involving Af^) and 
I is the unit operator. Furthermore we consider the angular-dependent term in this expansion 
to be a perturbation on the first term, Vo{r). For each value of J the Hamiltonian Hj will give 
a manifold of states which is the direct product of a manifold corresponding to various numbers 
of localized phonons being excited with the manifold of (2 J + 1) states having different values 
of mj. An important simplification is that spherical harmonics with I > 2J have no nonzero 
matrix elements in the manifold of states of angular momentum J. 

Note that apart from the kinetic energy, this Hamiltonian is a strictly local operator. Thus 
we solve the eigenvalue problem on a discrete mesh of points on a cube centered at the octahedral 
site when the wave function is required to vanish on the boundary of the cube. Each edge of the 
cube is taken to be [— L, L] with mesh point spacing of dL. In this scheme the wave function at 
each mesh point is a (2 J + l)-component vector. We are mainly concerned with the manifold 
J — and J = 1, in which case the problem is numerically not significantly harder than for a 
scalar problem. Even though the resulting matrix size is very large, it is a block band matrix 
and is very sparse. The numerical results reported here were obtained from L=1.65 A and 
dL = 0.075 A, which requires diagonalization of a matrix n x n where n — 273375). However 
we confirmed that a coarse mesh points with L — 1.2 A and = 0.17 A, gives almost the same 
results (where the matrix size is n — 14,739). The large sparse matrix eigenvalue problem is 
solved using the package ARPACK.-'^''' 

Since numerical results sometimes do not provide complete insight into the nature of the 
solutions, we have also used perturbation theory to understand the results. In this approach 
we treat Vj{r) in Eq. (3) as the perturbation. The unperturbed problem, apart from the 
additive energy BJ{J + 1) is thus that for translations of the spherical (J = 0) molecule. This 
spectrum is not too different from that of a three-dimensional harmonic oscillator. Accordingly, 
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to qualitatively interpret our more accurate numerical results, we apply perturbation theory 
in which we develop an effective Hamiltonian^® to describe the splitting of this manifold which 
is characterized by a value of J and of N, the total number of phonon excitations. (The 
perturbative effects due to coupling between manifolds of different J is negligibly small for 
hydrogen in Cgo-^) This effective Hamiltonian is a matrix of dimensionality D, where D = 
(2 J + 1)(A^ + 1){N + 2)/2 and schematically is of the form 

n{N, J) = BJ{J + 1)1 + Tiphonon + KT'^'^ - Vp^^^VP^^ , (4) 

where Tiphonon gives the energy of the various states with a total of N phonons. These energies 
are just those calculated for a (J = 0) molecule. Also Vj*^^*^ is the part of Vj which is diagonal 
with respect to the number of phonons, VP^^ is the part of Vj which is off-diagonal with respect 
to the number of phonons, and £ is the change in phonon energy caused by VP^^. 
This effective Hamiltonian is defined by its matrix elements as 

{N, a; J, M\n{J, N)\N, a'; J, M') = [BJ{J + 1) + En,^]5^^^,5m,m' 

+ j2{N,a\A^r'''i^)\N,a'){JM\Y,f~^''{n)\JM') 
1=1 

N'jtNl,l'=l fJ' 13=1 

{JM\Y,f-^mjf,){Jf,\Y,^,r^''{n)\JM') , (5) 

where kN — {N + 1){N + 2)/2 and the states with N phonons are labeled N, a, where a runs 
from 1 to k^. 

We now briefly discuss how the A^''s of Eq. (2) are obtained from the atom-atom potential 
between each H atom and each carbon atom. Here we will assume that the atom-atom potential 
is of the form F(|rj — rjj|), where and are the displacements of the ith carbon and of the 
H atom, respectively, relative to the center of the H2 molecule. For this form of potential, we 
show in Appendix A that 

A^ = 2nY: YrihT f (3a:^ - l)F {[r] - {p" + pxr,]"') dx , (6) 

where p is the separation between H atoms in the H2 molecule and the sum over i is over all 
relevant neighboring carbon atoms. It is instructive to expand this expression in p/r^, which 
yields 
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Af = Y.{f" - ^)yr(ii) + o(p'lrt), (7) 

where F" and F' are the second and first derivatives of F{ri). 

This expansion is good enough to reproduce most of the results discussed in this paper. 
We note that the A^'s (i.e. the orientational potential) are zero for a harmonic potential (i.e. 
F{ri) = \krf) because the prefactor [F" — F' /vj^ is zero. This can be also seen easily as follows. 
Assuming an atom-atom potential between each H atom and the C atoms in the adjacent Ceo 
molecules, we may write the potential of an H2 molecule as 

V{v- n) = K(r + Ipn) + K(r - Ifm) , (8) 

where is the potential of a single atom due to the entire octahedral cage in which it is 
confined, n is a unit vector along the axis of the molecule, and p is the separation between 
atoms in the molecule. As we shall see, the total potential is nearly isotropic. So we write 

Va{r)^lkr^ + Sr\ (9) 

When we substitute this into Eq. (8), we obtain the result 

V{r; n) = Voir) + 25[r^p^ cos^ - i] , (10) 

where 9r,n is the angle between the vectors r and n and Vq is independent of „. The point is 
that the orientationally dependent part of the interaction depends on the anharmonicity: for 
a purely harmonic and isotropic interaction T4, the total potential energy is independent of 
the molecular orientation. Thus we expect rotation-translation coupling to be weak. On the 
other hand in nanotubes, where the quadratic term is anisotropic, this coupling will be more 
important. 

III. ENERGY SPECTRUM OF A ( J = 0) H2 MOLECULE 

We start by considering the eigenvalue spectrum of Ho in which the orientational dependence 
of the potential is neglected. In this approximation, apart from the additive constant BJ{J+1), 
the total energy (rotational plus translational) is the same as that of a ( J = 0) molecule. For 
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most purposes a (J = 0) molecule may be considered to be a spherical molecule because 
the orientational wave function ¥^(0,) is uniform over all orientations. Each eigenfunction of 
Hoi J) is the product of a rotational function taken from the manifold of 2 J + 1 degenerate 
orientational wave functions and a translational wave function which represents an eigenfunction 
for a spherical molecule confined to a cage. These translational wave functions satisfy 



where Vo{r) is the potential discussed in Appendix A. The index k labels states which we 
might otherwise label by three indices, each quantum number characterizing the number of 
excitations in each direction. Note that these unperturbed solutions do not involve the coupling 
between rotations and translations. As discussed above, these eigenfunctions were obtained by 
converting the continuum equation (11) into a into a discrete equation on a mesh of points and 
solving the resulting matrix eigenvalue problem using a sparse matrix routine. 

Since it happens that the energy levels and eigenfunctions we obtained numerically are not 
qualitatively different from those of a spherical harmonic oscillator, we first study the energy 
spectrum as perturbations, S, k, and A, are sequentially turned on in the following potential: 



Fig. 1 shows the evolution of the energy spectrum as perturbations are sequentially introduced 
which take the spherical harmonic oscillator into the actual lower symmetry of a molecule in 
an octahedral interstitial site. In the left-most panel we show the energy levels for a spherical 
harmonic oscillator, with hcu adjusted to correspond to the single-phonon levels of H2 in an 
octahedral interstitial site in Ceo- Note that the levels are highly degenerate because the energy 
depends only on the total number of quanta. The symmetry of the Hamiltonian is ^7(3), the 
group of unitary three dimensional matrices. We now add to this potential an anharmonic term 
of the form Sr^. This perturbation lowers the symmetry to that of the rotation group in three 
dimensions. As is well known, each eigenfunction in a generic spherically symmetric potential 
can be labeled by the magnitude of the orbital angular momentum, L. Thus the single phonon 
levels are unsplit by this anharmonic perturbation and are now labeled as angular momentum 
L = 1 states, whereas the two phonon levels split into a manifold of five L = 2 states and one 




(11) 




|r^) + X{xy + yz + zx). 



(12) 
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L = state and similarly for states with more than two phonons. In Fig. 1 we have taken the 
constant S to be that which best describes the anharmonicity of H2 in Cgo- The energies of the 
perturbed levels are given in Table I. 

Next, we consider what happens when the spherical oscillator potential is augmented by a 
cubic symmetry potential of the form /t(a;*^ + + z*^ — This potential is appropriate for a 

spherical molecule in an octahedral interstitial when the Ceo are orientationally disordered and 
have an Fm3m crystal structure/ The degeneracy associated with spherical symmetry is lifted, 
but as shown here one retains cubic symmetry, so the three one-phonon states which transform 
as X, y, and z are degenerate. The two-phonon states are of three different symmetries. One 
{t2g) transforms like xy, xz, and yz. This is the lowest level. The next highest level is the 
s-wave symmetric combination which transforms like x'^ + y"^ + z'^. Then one has a doublet 
of d-wave (e2g) symmetry. This classification scheme is continued in the higher-energy levels. 
Although we are not dealing with harmonic phonons, it is still useful to consider manifolds 
characterized by the quantum numbers J and A^, which are respectively the rotational angular 
momentum and the total number of phonons, at least up to = 3. Quantitative results are 
given in Table II. 

Finally, in the rightmost panel of Fig. 1 we show the further reduction in degeneracy 
which occurs when the octahedral interstitial is surrounded by Ceo molecules which have the 
long range order associated with the Pa3 crystal structure. ^-^^ In this case, each interstitial is 
uniaxial (with symmetry Se) rather than octahedral. Accordingly, we introduce a potential of 
the form X{xy + yz + zx) = |A(3^^ — r^), where the ^ axis is taken to lie along the three-fold 
axis of the interstitial site. There are four symmetry related interstitial sites, each of which 
has its three-fold axis along a different [1, 1, 1] direction. The resolution of degeneracy in the 
presence of this uniaxial perturbation is also given in Table II. In all these cases, no interactions 
between rotations and translations are involved. 

We have solved the eigenvalue problem of Eq. (11) on a mesh of points and obtained 
the results given in Table III. Results labeled "Octahedral" are those for the orientationally 
disordered phase, where each Ceo molecule is replaced by a sphere of carbon atoms as is discussed 
in Ref. 5. Since these numerical results lead to manifolds of energy levels associated with a given 
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number of phonons and the degeneracies of these manifolds are as expected from our general 
discussion above, we conclude that the potential seen by a spherical H2 molecule in the low- 
lying phonon levels is not very different from that of a spherical harmonic oscillator. However 
as noted in Ref. 5, the effective harmonic potential must be taken to be a self-consistently 
renormalized potential to take account of the larger zero-point motion. 



IV. ENERGY SPECTRUM OF A J = 1 MOLECULE 

We now discuss the energy spectrum of an ortho molecule with {J — 1). As we have seen for 
{J — 0) molecules, our numerical results indicate that for N up to, say, three, one can clearly 
identify the manifold of N phonons. We therefore discuss the systematics of these manifolds. 

A. Zero-Phonon Manifold 

We first consider the case of J = 1 with = phonons. This manifold is described by the 
effective Hamiltonian 

H=[2B + AE]I + S[J^, - |] . (13) 

The splitting 6 must be zero when Ceo is orient at ionally disordered. From Eq. (5) one sees 
that because the spherical harmonics are traceless, the average energy shift, AE, has nonzero 
contributions only from terms which involve coupling to excited phonon states. (In Ref. 5 a 
neghgibly small shift was found due to effects off-diagonal in J which we ignore here.) Prom 
Eq. (5) we find that 

^E^-l E ^^:«El(0,l|^5(r)|A^,a)n(l(/. + r)|y,^(Q)|l/.)r . (14) 

To implement this equation, we first construct AY^{r) as discussed in Eq. (6). Then matrix 
elements of A^{v) are taken between phonon states for a J = molecule which we obtained 
previously and which are labeled N,a (|0, 1) being the phonon ground state). Thereby we 
obtained the results given in Tables IV and V. In Eq. (14) the matrix elements of spherical 
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harmonics Y^{fl) are taken between orientational states labeled by J and J^. To evaluate AE 
we use 

Y.mH + r)\Y,-mM\'-^ (15) 

so that 

^E = -1KI E ^i^,aEK0'l|^2(r)|iV,«)r (16) 

In appendix B we give a model calculation of an H2 molecule in a spherical cavity from 
which we evaluate Eq. (16) to give AE — — 0.14meV. In this calculation the translational 
wave functions are assumed to be those of a harmonic oscillator with (r^) = 0.1875^^. As 
noted, the result is very sensitive to the value used for (r^). For octahedral symmetry (i. e. 
for orientationally disordered Cqo) we evaluate Eq. (16) using the data in Table IV. Thereby 
we find a shift AE = —0.133 meV. The same approach using our numerical solutions for the 
phonon states of a (J = 0) molecule for the orientationally ordered Pa3 phase yields the result, 
AE — —0.141 meV , compared to the experimental value^ AE — —0.35 meV. Again we mention 
that a small change in parameters could easily lead to a much larger calculated value of AE. 
Prom the numerical solution for the three component wave function of a ( J = 1) molecule on a 
mesh of points, we obtained the value AE = —0.16 for the Pa3 phase. The various numerical 
results for AE are summarized in Table VI. 

Prom Eq. (5) we also find the sphtting (in the Pa3 phase) to be 

3(0,1|A0|0,1) 



+J_VE(e)-i 
207r4" ^ ^ 



1(0, l\Al{r)\e)\' + 1(0, l\Al{r)\e)\' - 2|(0, l|^^(r)|e) 



(17) 



where the quantization axis is taken to lie along the threefold axis of symmetry of the interstitial 
site. Using the matrix elements given in Table IV, we find that the contribution to the splitting 
S comes almost exclusively from the diagonal term (0, 1| A2IO, 1) and we obtain the results listed 
in Table VI. 
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B. One-Phonon Manifold 

1. Numerical Results 

Next we consider the manifold J = 1 with = 1 phonon. Again only YJ*" with / = 2 
contributes, so that we may write 

l)a^;aV' = {Mn{N = 1, J = l)|aV) 

= [2B + Eo + /lu;.] W^a,a' + {la\A^-'''ir)\la')MYr'''imifJ'') 
-EE E^^ ^{la\A^,-'\r)\NW'){N'a''\A{-^\r)\la'^ 

x(i/.|rr^"(i])|i/)(i/|rr-^'(i])|V) . (18) 

This is a 9 X 9 matrix, which gives the 9 ( J = 1, A?" = 1) levels. Since the first term gives 
rise to the removal of degeneracy expected from group theoretical considerations, we did not 
include the second term in our numerical evaluations. This procedure was sufficiently accurate 
to provide a useful check on the validity of the more accurate numerical solutions for the three- 
component wave functions of our set of mesh points. In Table VII these numerical results ( "full 
mesh") are given and are compared to the results using perturbation theory, as in Eq. (18). 
As can be seen, the two approaches yield quite compatible results. 

2. Qualitative Remarks 

Some additional comments on Eq. (18) are in order. The first fine of this equation gives 
the energy at first order in perturbation theory. At this order the wave function remains a 
product of the spatial ground-state spatial wave function for a J = molecule times a J = 1 
rotational wave function. At this level of approximation there is no dynamical coupling between 
translation and rotation. In second order perturbation theory we see that admixtures of two 
phonon states which are multiplied by different rotational states are introduced. For example, 
consider the situation when the molecule is in a uniaxial symmetry site and let \X), \Y), and 
\Z) be the J — 1 states for which respectively Jx, Jy, and are zero. If, for simplicity, we 
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assume that the unperturbed spatial wave function is spherically symmetric, then the state 
which without perturbation was 

Co|Z)e-'''/(^"') (19) 
where Cq is a normalization constant, is now 

where Ci ~ Co and C2 is small compared to Ci. The point is that this formulation allows the 
molecule to change its orientational state as it translates. For H2 in Ceo this effect is small, 
however, in less symmetrical cavities, as in nanotubes^^, this effect can become more important. 

3. Group Theoretical Analysis 

In Fig. 2 we show the influence of roton-phonon coupling and local site symmetry on the 
energy spectrum of the one phonon-(J = 1) manifold. At the far left we start from the case 
of highest symmetry when the phonon and rotations separately have complete rotational in- 
variance and no phonon-roton coupling is present. In this case the manifold of 9 states (3 
one-phonon states (g) 3 (J = 1) states) is completely degenerate. When roton-phonon cou- 
pling is included (but the environment is still spherically symmetric) we have overall rotational 
invariance and the resulting cigenstates are characterized by their total angular momentum 
K. The roton-phonon coupling causes states with different K to have different energy. (The 
size of the splitting shown in the figure is adjusted to agree with the center of gravity of the 
appropriate levels for cubic site symmetry.) 

The two right-hand columns pertain to the situation when a (J = 1) hydrogen molecule 
occupies the octahedral interstitial site of Cqo- When the Ceo molecules are orientationally 
disordered the interstitial site has Oh symmetry and we consider that case first. Use of the 
character tables for the Oh group indicates that the original 9 dimensional reducible represen- 
tation r is decomposed into irreducible representation of the Oh group as 

r = Tag e Tig ®Eg®A^. (21) 



14 



and the basis functions associated with these irreducible representations are given in Table 
VIII. As mentioned above, for temperatures below about 260K, the Cgo molecules order into 
a structure of crystal symmetry Pa3^^^ in which case the formerly octahedral interstitial has 
the lower Se symmetry. Use of the relevant character table shows that now 

r = 3Ag e 3Eg e 3^;* , (22) 

where Eg is a complex one dimensional representation and E* is its complex conjugate partner. 
The basis functions associated with these irreducible representations are given in Table VIII. 
The most important conclusion from this analysis is that the energy eigenfunctions are not 
simply products of translational and rotation wave functions, but instead are linear combina- 
tions of such products. This type of wavefunction reflects the fact that symmetry operations 
act simultaneously on the position and the orientation of a molecule. 

To emphasize this fact we give, in Fig. 3, a pictorial representations of the translation- 
rotation wavef unctions. This representation is to be interpreted as follows. We know that 
the rotational wave functions for a free J = 1 molecule can be taken to be analogs of Px, Py 
and Pz functions and we will label these rotational wavefunctions as X, Y, or Z. For instance 
\X) ~ sin ^ cos 0, \Y) ~ sin ^ sin 0, and \Z) ~ cos^. These wave functions have two lobes, one 
positive the other negative aligned along the axis associated with their state label. When such 
a rotational function is multiplied by a one phonon function in the a direction, (|a;) denotes a 
wave function for a single phonon excitation in the x-direction), the total wave function will 
be an odd function of a. Thus the wave function \xX) is an odd function of x and is therefore 
depicted by two Px functions, one at positive x and another at —x with the signs of the two lobes 
changed. For simplicity in the figures we show only those functions which have appropriate 
dependence in the plane of the paper, which is taken to be the x-y plane. At the upper left 
we show an xy-like function. It has two other t2g symmetry partners which are xz-\ike and 
yz-like. At the upper right we show one of the two Eg functions which is x^ — y^-like. These 
five t2g and Eg functions comprise the manifold of total angular momentum K — 2 states. 
Within spherical symmetry all five of these states are degenerate in energy. In the lower left 
of Fig. 3 we show the 2;-like function of tig symmetry. Its other two partners are obtained by 
cyclically permuting a;, y, and z. These three functions comprise the manifold of total angular 
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momentum K = 1 states, which transform under rotation as a vector. Finally at the lower 
right we show the angular momentum K — state. Thus in spherical symmetry, the nine J — 1 
single phonon states give rise to three distinct energy manifolds which have degeneracies 1, 3, 
and 5, corresponding respectively to total angular momentum K — 0, K — 1, and K — 2. 

The simplest classical arguments do not reproduce the above results. For instance, one 
might argue that translation can occur equivalently along either of three equivalent coordinate 
axes. In each case, one can have the molecule oriented along the axis of translational motion 
or perpendicular to that axis. This argument would suggest that the nine levels break into a 
three-fold degenerate energy level in which the molecule is oriented longitudinally and a six-fold 
level in which the molecule is oriented transversely. This discussion shows that it is essential 
to treat the translation-rotation problem quantum mechanically to get the correct degeneracy. 

V. NEUTRON SPECTRUM 

In the experimental study of Fitzgerald et. al.^, neutron energy loss spectrum of H2 trapped 
in Ceo was measured with an energy resolution of 0.3 meV. The spectrum shows surprisingly 
rich features. However the origin of these features were not successfully identified in detail. 
Since the observed neutron spectrum is a direct probe of the intermolecular potential between 
H2 molecules and Cgo host lattice, it is very important to see if available atom-atom potentials 
can give a spectrum which is similar to the experimental data. A suitable analysis of the 
high resolution inelastic neutron scattering data in Ref. 5 should, in principle, give a detailed 
information about the intermolecular potential between H2 molecules and the host lattice. 

Figure 4 illustrates several possible transitions, involving both rotational and vibrational 
excitations, that could be observed in a neutron scattering experiment for H2 in solid Ceo- In 
order to estimate the intensities of these transitions and the corresponding neutron spectrum, in 
Appendix C we derive the inelastic neutron cross section for trapped H2 molecules in a powder 
sample at low temperature. Below we discuss the contribution to the total neutron spectrum 
from each of these transitions, labeled as Ta, Te, and then compare the calculated spectrum 
with experimental data using various atom-atom potentials. 

We start with the transitions involving phonon creation in para hydrogen, as shown by 
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Td in Fig. 4. Because of the spin-dependent interaction between the proton and the neutron, 
processes in which a para molecule is not converted to an ortho molecule are forbidden, or 
more correctly speaking, are proportional to the coherent cross section b, which is very small 
compared to incoherent cross section b'. Hence, the transition To will not have a noticeable 
contribution to the total cross section. 

We next discuss the contribution to the total spectrum from processes in which either a 
(J = 0) molecule is converted to (J = 1) molecule (para-ortho conversion, labeled as Ta in 
Fig. 4) or a single phonon is created (as shown by Tb). Our calculations presented so far 
indicate that both processes will give features around 14 meV in neutron spectrum. 

In the Appendix C we find that the cross section due to para-ortho conversion, (indicated 
by the subscript — > 1 is given by 

^ = fiv|(l-x)[6'j,(i«:p)]V^^WX:5[^L-(i^c + i^..)] , (23) 

where is the total number of H2 molecules, k (k') is wave vector of the incident (scattered) 

neutron, k. = k' — k, a; is the fraction of H2 molecules which are ortho (odd J) molecules, p is 

the separation between protons in the H2 molecule, b' is the spin-dependent cross section in the 

proton-neutron pseudopotential, jn is the nth order spherical Bessel function, and W{k) is the 

Debye- Waller factor which we take to be ^k^{u^). Also, is the energy loss of the neutron, 

and Ec + Em is the para-ortho conversion energy when the final state of the ortho has = m. 

Similarly, the cross section due to ortho-para conversion, q^q^ ) has the same expression 

/ 1^0 

as Eq. (23) but now the factor {1 — x) is replaced by x. Hence the ratio of the total cross section 
for ortho to para conversion to that of para to ortho conversion is (1 — x) to x, where x is the 
ortho concentration. Normally the ratio of energy gain to energy loss cross sections follows the 
Boltzmann factor. Here, the populations are set by x rather than by the temperature. 

We now discuss the cross section due to phonon creation on a (J = 1) molecule (indicated 
here by the subscript 1^1), which is calculated in the Appendix C. These transitions are 
shown as Tb in Fig. 4. The result requires a knowledge of the translation-rotation wave function 
of the H2 molecule. We find that 

^)^^_ = iv4'^i:/lV. (24) 
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where the cross section S'^i^^ „ are given in Eqs. (C34), (C35), and (C43) of Appendix C. 

In Fig. 5 Tb represents the transitions from the ( J = 1, = 0) levels to the manifold of 
nine energy levels of(J=l,A" = l). Accordingly, we expect several transitions with nonzero 
amplitude and thus rich features in the total neutron cross section. 

Also one may consider the cross section integrated over energy, which is a usefull quantity 
to indicate the relative strength of the different transitions discussed above. The ratio, rp, of 
the integrated neutron energy-loss cross section for phonon creation (process Tb in Fig. 4) 
divided by that for para-ortho conversion (process Ta in Fig. 4) was found in the Appendix C 
to be 

This ratio is plotted as a function of k for x = 3/4 in Fig. 6. Since this ratio is of order unity, 
the energy loss spectrum will display features due to both phonons and para-ortho conversion. 

In appendix C we also calculate the zero phonon ortho cross section for the transition shown 
as Te in Fig. 4. Since the (J = 1) levels are in thermal equilibrium, in this case the ratio of the 
cross sections at energy gain to those of energy loss does satisfy detailed balance. The ratio of 
the total cross section (counting both energy loss and energy gain) for transitions within the 
{J — 1) ground manifold to that due to para to ortho conversion was found at zero temperature 
to be 

rj=i{T = 0) = — — ^ — - . 26) 

1-x 15ji(2«;p) 

Figure 6 shows that this ratio is quite small and therefore experimental observation of this 
transition (i.e. Te) would be very difficult. 

Figure 7 shows the neutron energy loss spectrum and the calculated total spectrum using 
the same potential, so called WS77 model, used by FitzGerald et. af. Even though the 
calculated spectum is wider than the experimental spectrum, it is still possible to make a one 
to one correspondence between calculation and experiment as is shown by arrows in Fig. 7. 
The top curve in this figure shows what the spectrum looks like if the orientational part of the 
potential is scaled by about half. The agreement between the data and calculations is somewhat 
better after this arbitrary scaling, indicating that the potential used is too anisotropic for the 
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center of mass motion of H2 molecule but at the same time it is too weak for the orientational 
part of H2 (because it gives too small a result for 5, the sphtting of the {J — 1) ground 
manifold) . 

We also tested other potentials commonly used in the literature and these results are shown 
in Fig. 8. The top curve is from Novaco's 6-12 potential which was developed to study hydrogen 
on graphite."^ Clearly this potential gives too low phonon energies and too little splitting of the 
(J = 1) levels for H2 in sohd Ceo- The other two curves in Fig. 8 are 6-exp potentials tabulated 
in Ref. 15. The spectrum from these potentials does not agree with experiment either. We 
also searched the potential parameters A and B for 6-12 and A, B, and C for 6-exp types of 
potentials. However wc were not able to improve the fit to experiment using these atom-atom 
potentials. Hence, it seems that simple atom-atom potential does not describe the details of 
the H2-C60 interaction well. It is an open and important question to find a better potential 
can reproduce the experimental spectrum better. It is also of important to see how well the 
potentials obtained from density functional theory within local density approximation will do. 

Finally, in addition to features observed near 14 meV, Fitzgerald et al.^ also observed a 
feature at about 28 meV in the energy gain time-of-flight spectrum. This energy is about twice 
that of either the para-ortho conversion energy or the energy of the translational phonon for 
an H2 molecule in an octahedral interstitial site. Clearly this feature represents the energy of 
two excitations, but it was not clear whether these would be two phonons, one phonon and 
one ortho-para transition or two ortho-para transitions. In Fig. 9 wc show the temperature 
dependence of the total intensity of this feature. This temperature dependence follows a thermal 
activation with an energy of about 14 meV. Thus the initial state must consist of one thermally 
excited phonon and the transition observed destroys one thermal phonon and converts one 
ortho molecule (which occurs with temperature-independent probability x) to a para molecule, 
thus giving the observed energy of about 28 meV. 

VI. CONCLUSIONS 

Our results indicate that the coupled phonon-roton problem is a rich one. For the light 
molecules of H2 where the sphtting of rotational levels is large compared to most lattice vibra- 
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tional modes, one is in a so-called weak-coupling limit where the interactions between rotons 
and translational phonons can be treated perturbatively. Even at this level one recovers an 
interesting structure. 

Needless to say, we hope that the calculations in this paper will motivate more detailed 
experiments at higher resolution to elucidate the structure of the roton-phonon spectrum. 
We may summarize our main conclusions as follows: 

• We have presented a systematic perturbative approach to the calculation of the roton- 
phonon spectrum of hydrogen molecules in confined geometry. Our calculations agree with the 
group theoretical analysis for the geometries considered here. 

• In a general way, the techniques of this paper (use of the atom-atom potential combined 
with perturbation theory) may prove useful to treat hydrogen molecules in other confined 
geometries, in particular in or on nanotubes. We are currently analyzing this situation. 

• We give the first calculation of the expected energy loss spectrum from hydrogen in Ceo 
in the energy range where phonons and para-ortho conversion both are important. We find 
that none of the traditional 6-12 and 6-exp types of potentials give good results for the detailed 
energy dependence of the observed phonon spectrum, although the WS77 potential^^ we used 
was definitely the most satisfactory. It is a theoretical challenge to determine a potential which 
fully reproduces the observed spectrum. 

• We identify the feature at 28 meV in the energy gain spectrum as consisting of conversion 
of one ortho molecule to a para molecule combined with annihilation of a single phonon. This 
identification is uniquely indicated by the temperature dependence of this feature. 
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APPENDIX A: ATOM- ATOM POTENTIAL AND AfJ 



We model the interactions between the hydrogen molecule and the surrounding cage of Cgo 
molecules using an atom-atom potential. Unless otherwise indicated, all the results reported 
in this paper are obtained from the same potential —A/r^ + B exp{—Cr) that is used in Ref. 
5 (where A = 5.94 eV A^, B = 678.2 eV, and C = 3.67 A-^). 

We consider two cases depending on the orientational state of Cqo molecules. When the 
molecules in the surrounding cage are orientationally disordered, we distribute the carbon cen- 
ters uniformly over the surface of a sphere. This case corresponds to the octahedral symmetry 
discussed in the text. The resulting integration of the atom-atom potential over a spherical 
surface is done analytically in Ref. 5 and therefore is not given here. For the Pa3 symmetry, 
the Ceo molecules arc rotated according to their Pa3 settings and then the total potential and 
Af'^s are calculated on a mesh points of a cube centered at the octahedral site. Below we derive 
a convenient way to obtain the ^^'s from the atom-atom potential. 

We write the potential, Vh-c between a single H atom and a single C atom as 

Vh-c = F{r) , (Al) 

where r is the displacement of the C atom relative to to the H atom. Thus the interaction, 
y(H2) of a H2 molecule with a C atom can be written as 

V{R,) = Y: nir' + \p' + cjpr ■ , (A2) 

cr=±l 

where now r is the displacement of the C atom relative to the center of the H2 molecule whose 
atoms are at positions ±|pn, where n is a unit vector specifying the orientation of the molecular 
axis of the H2 molecule. Then 

= E E / y27inrF{[r^ + + apv, ■ n]'/')dn , (A3) 

a- i ■' 

where (Kt indicates an integration over all orientations of fi and the sum over i is over all 
neighboring carbon atoms. We use this to get 

AI^V,^2^Y: £ + 1/ + prixY'^)dx . (A4) 
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To get the Af^^s for / > we write 

E nir- + Ip' + <^pr^ ■ n]'^') = E B2L{r,)Y,',{er,n) , (A5) 

cr L 

where 9r,n is the angle between the vectors rj and n. We have that 

S2L(n) =27r$: r F{[r'^ + \p'' + apriCOseY/^)Y^^^{9)sm0d0 . (A6) 
Now substitute Eq. (A5) into Eq. (A3) to get 

= E / ^2?(^)* E 52L(rOy2°L(^.,n)rf^^ ■ (A7) 
i •' L 

Using the addition theorem for spherical harmonics^'^ we have 

= E / >^2?(^)*E^2L(r.)VV(2^+l)E^2LN>TL(i-^)*^^^ 



= E^2Kn)VV(2^ + l)1^27(r.)* ■ (A8) 

i 

For we get 

A- = 27r^yr(rO* /' (3x^ - + + npx]'^')dx . (A9) 

For power-law functions F (i.e. F ~ l/r^"), this integral can be done analytically. 

APPENDIX B: SPHERICAL CAVITY 

In this appendix we apply our formalism discussed in the text to a simple toy model to 
utilize the main physics of quantum roton-phonon dynamics of a H2 molecule confined in a 
spherically symmetric cavity for which the potential differs perturbatively from harmonic. 

1. Orientationally Dependent Potential 

For a diatomic molecule for which the spherical part of the potential is that of an isotropic 
spherical oscillator, we take the orientationally dependent part of the potential to be 

U{r,fl) = f{r) {x sinO cos(f> + y sinO sincf) + z cosOy — ^ , (Bl) 
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where i. = x/r, y = y/r, and z = z/r. 

This potential can be written in the canonical form of Eq. (2) with 



AT = %f{r)Yr{fr . (B2) 



2. Energy Shift of the J = 1 Manifold 

To evaluate Eq. (16) for the shift in the center of gravity of the J — 1 zero phonon levels, 
we need the wave function of the ground state, namely 

V'o(r) = ae-'-'/^^'^') . (B3) 

where a — a~^/'^{2n)~^^'^. Also we can write the two phonon excited states (for spherical 
symmetry) as 

- /3(r/a)Vr(r)e-^'/(^'^') , (B4) 



for m = —2, —1, 0, 1, 2 and for angular momentum L = 2 and where P = a'^^"^ y 2 /15{2tt)^^^^ . 
The sixth two-phonon state is an s-wave state, whose wave function we do not need. We will 
also need the L — 2 (d-wave) four phonon states which are 

4f='^ = liir/a)' - 7{r/ar]Yr{r)e-^'/^'^'^ , (B5) 

for TO = —2, —1, 0, 1, 2 and for angular momentum L — 2 and where 7 = (7~^/^-y/l/105(27r)~^/^. 
Now we assume the following polynomial fit for /(r): 

/(r) = 72(r/<7)2 ^ 74(r/<7)^ + ^e{r/af . (B6) 

We then find 

l(0|A™(r)|4fj'))| = 4^^172 + 774 + 63761 = (B7) 

and 

|(0|A-(r)|V^if=^))| = 4^174 + 18761 (B8) 
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Thus we have 

_ 4 (72 + 774 + 6376)^ 28 (74 + 1876)^ 

15 hu 15 hu ' ^ ' 

Prom numerics we get 

Al(z) = 0.07226r^ + 0.0348r^ . (BIO) 
For convenience we take a — 0.25A. Then A2 is of the form of Eq. (B6) with 

74 = 0.07226(7^ , 76 = 0.0348(7^ . (Ell) 

74 = 0.07226/256 = 0.282 meV (B12) 
76 = 0.0348/4096 = 0.0084 meV . (B13) 

we evaluate the first and second terms on the right-hand side of Eq. (B9) to be 0.115 meV and 
0.025 meV, respectively. Thus we have 

AE = 0.14meV . (B14) 

Notice that AE is a very strong function of a. For instance, if you take a — 0.26^4, you get 
74 = 0.33 meV and 76 = 0.011 meV in which case AE — 0.19 meV. 

3. Effect of Translation-Rotation Coupling on the J = 1 One Phonon Manifold 

In the absence of coupling between translations and rotations we characterize the single- 
phonon states by phonon angular momentum Lp, so that 

\Lp = ±1)^T (^^) ^o(r) , \Lp = 0) = {z/a)M^) . (B15) 

We now wish to include the effect of the perturbation of the form of Eq. (2) when the coefficients 
are as in Eq. (B2), with /(r) given by Eq. (B6). We know that states are now characterized 
by the total angular momentum K = Lp + J. So the energy of the K — 2 manifold is given by 
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E{K = 2) = {Lp = 1, J, = l|f/(r, n)\Lp = 1, J, = 1) 



1 Jdr 


x—iy 
a 




15 


Jdr 


x—iy 
a 





;j, = i|(3j,2-2)|j, = i) 



2/drr2/(r)e-'^'/(2^') 



(B16) 



225(72 /dre-'^V{2a2) ' 
where /(r) and F are defined in Eq. (B6) and Eq. (B7), respectively. 

Similarly, one can evaluate the Hamiltonian in the manifold of states 0i = \Lp = 1, = 0) 
and (f)2 — \Lp — 0, Jz — 1). The Hamiltonian matrix in this basis is found to be 



J_ 2 
15 5 



15 



(B17) 



This matrix has an eigenvalue which is associated with the K — 2, — I state and the 



3^ • 



new eigenvalue for the K = 1 manifold, E{K = 1) : 

Similarly, one can evaluate the Hamiltonian in the manifold of states (pi = \Lp = 1, Jz = 
— 1), 4>2 — \Lp = 0, Jz = 0) and 03 — \Lp — —1, = 1). The Hamiltonian matrix in this basis 
is found to be 



n 



2_ 
15 



2 4 
'5 5 



2 _8_ 

5 15 



2 A 
'5 15 



(B18) 



In this manifold we reproduce the eigenvalues for X = 2 and K — 1. The new eigenvalue is 
E{K = 0) = |F. 

APPENDIX C: NEUTRON SCATTERING CROSS SECTION 
1. GENERAL FORMULATION 



Following Elliott and Hartmann^^, we write 



mo 



-5(r-r„) 



b + b'(a ■ I) 



(CI) 



where mo is the neutron mass, cr is the neutron spin, r is the coordinate of the proton, r„ the 
coordinate of the neutron, / is the proton spin, and b and b' are the coherent and incoherent 



25 



scattering lengths for this scattering. Since b is very small, we drop that term from now on. 
The differential scattering cross section is 

k' 



dQdE k ^ 



^PiS{E-Ei + Ef)\V\^ , 



(C2) 



if 



where E = n'[{k'Y -k^]/(2mo) and the sum is over all states of the system. Here the potential 
is 



(C3) 



where k. = k' — k, j labels molecules and a = 1,2 the atoms within a molecule. Performing the 
sum over a we get 



(C4) 



where Ij = 1^1+1^2. The first term acts only on ortho molecules because for paras the total spin 
is zero and the second term causes transitions between ortho and para molecules. So we write 
the scattering cross section as the sum of three terms, the first of which represents scattering 
from an ortho molecule and others ortho-para conversion or the reverse. Thus 

k' 



NxSi^i + NxSi^o + N{1 - x)So^i 



(C5) 



dfldE k 

where N is the total number of molecules and x is the ortho concentration. Because there are 
no correlations between nuclear spins each cross section is actually a sum over cross sections 
for each molecule: 



where j labels the molecule, /3 is — > 1, etc., and 



(C6) 



So^lj = J2 + Ef) 

Ji=0,Jf=l 

Si^oj = E Pi^iE -Ei + Ef) 

Ji=l,J^=0 

Si-^i,J= E mE-E, + Ef) 

Ji=i,jf=i 



(/|6V'^-^^sin(i^.-p)[cT-(I,,-I,2)]|^)T 
(/|6V'^■^^sin(iK•p)[o-•(I,l-I,2)]|^)T 



(/|6V^-^^cos(iK-p)[<T.I]|i)T 



(C7) 
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where the sums are over states for the fixed species (ortho or para) of molecule as indicated 
and the subscript T indicates that the wavef unctions include nuclear spin functions. 

Now we perform the sum over the spin states of the neutron and proton to obtain the results 

{f\e"^-^^smC^K-p)\i{ 



Ji=0,Jf=l 



Ji=l,Jf=0 



Si- 



>hj - 2 



E mE-E,+Ej] 



{f\e^'^-^^smaK-p)\i) 



(/|e^'*-^^cos(i/*-p)|i) 



Ji=i,jf=i 

where now states |/) and \i) no longer include nuclear spin wavefunctions. We write 



(0) 



e e 



(C8) 



(C9) 



where R^'^-' is the equilibrium value of R-,- and W ~ |k^(-u^) = {u^ + Uy + u'j.) . Since spherical 
harmonics of degree higher than two do not affect the manifolds of J = or J = 1, we now use 



cos(|^c • p) = joC^Kp) - 47rj2C^Kp)^Y,^{krY,"{p) 



(CIO) 



and 



sin(i« • p) = AnhaKp) E Y,^{krY;\p 



(Cll) 



We expand in displacements to get 

So^^j ^ {AttYA Y: PiS{E-E, + Ef){f\[l + iK.nj]j:Y,^{krY,^{p)\i) 

Ji=0,Jf=l v 

x(i| n,]Y.Yr{k)Yapr\f) , (C12) 



Ji=\,Jf=0 

En'^(^)*n^(p)K)(^| [1 - • n,]Y.Yr{k)Yrmf) , 



(C13) 



where A = le-^^[h' ji{\Kp)f and B = \e-''^[h' n{\Kp)f and 
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Ji=i,jf=i 



X(/| [1 + i/C-Uj] 

x(i| [1 — iK. • Uj] 



I/), 



(C14) 



Since the phonon energy is much larger than the orientational energy, we may classify 
contributions according the number of phonons that are involved. In the notation of Eq. (C6) 
we write 



(C15) 



where S^'^ corresponds to a zero-phonon process and S^^ to a process in which one phonon is 
created or destroyed. Thus 

Ji=0,Jf=l tiv 



4^i-(47r)2yl E Pi5{E-E, + Ef) 



Ji=0,Jf=l 

X E(/i<n'^(p)K)(^i«;n"(p)i/)'^««:/3n'^(^)*n"(A) . 



(C17) 



Here and below we use spherical components of a vector v: ^±1 = =F(f^a; =t iVy) / \/2 and vq = Vz- 
Similarly 



sf\ = {^^fB E mE -E,+ Ef) Y^ifiyrimmmmi^yya^) , (cis) 



Ji=l,Jf=0 



HI/ 



Si% = {iTTYB E Pi^iE -E, + Ef) 



Ji=l,Jf=0 

X E(/i<n'^(p)K)(^i«;n"(p)i/)'^"'^/3n'^('^)*n"(^) , 



(C19) 



5Si = (47r)^c E P.5(£;-£;, + £;^)E(/l>T(p)l^)(i|>T(p)l/m'^(^)%'^(^)*, (C20) 



ji=i,jf=i 



where C = ie-2^[6'i2(|Kp)]^ and 
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J,=l,.Jf=l a/3 
Ji=l,Jf=l //a/3 

+A7rD, J2 -E, + Ef) Y^{f\ul\i){i\u*pY,^(p)\f)KaKpY,^(k)* 

,Ji=l,J f=l fiaB 

+{A7:fD, Pr^iE -E, + Ef) ^ (/|<r2'^(p)*N)(i|4ir(p)|/)«:a«/3i^2"('i)*i"2'^('^) 

Ji=l,Jf=l nvaf) 

= E ^Sl;n , (C21) 
n=l 

where — \{h')'^e~'^^ jQ{^Kp)'^~'^j2{\np)^ and S^l^-^.^ is the contribution to the one-phonon 
ortho cross section from the nth term in the first equahty. 

Note the existence of terms in which a phonon and a roton are created, the system evolves 
and finally a phonon is destroyed. This type of process can only occur when the system 
supports roton-phonon interactions. All the terms <Si!^i;„ correspond approximately to the 
phonon energy. 

2. POWDER AVERAGE AT LOW TEMPERATURE 

Here we restrict attention to the energy loss spectrum at low temperature when there are no 
thermal phonons present. Also, we now take the powder average. This corresponds to actual 
experimental situation in Ref. 5, but would also be a reasonable approximation to take account 
of the differently oriented symmetry axes of the various octahedral interstitial sites. Below we 
calculate the cross sections for the following processes; (a) energy loss by conversion, (b) energy 
gain by conversion, (c) single phonon energy loss, and finally (d) zero phonon transition from 
( J = 1, M = 0) to ( J = 1, M = ±1). 

a. Energy Loss by Conversion 

We have 

{SZ,) = (47r)^ Y Pi^iE -E, + Ef) YifmpmmTpTlf) , (C22) 

Ji=0,Jf=l tJ. 
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where ( ) indicates a powder average. The initial state is the J = 0, = zero-phonon state, 
which we write as 



V^i = ^Q(r)|r;J = 0;J, = 0) , 



(C23) 



where Cj(r) is the amphtude of the wave function at the mesh point r. The final state is a J = 1 
zero-phonon state, which we similarly write as 



^f,m = E c/,m(r, M)|r; J = 1; J, = M) 



(C24) 



r,M 



and whose energy relative to the J — state is 
If El = —E is the energy loss, we may write 

(sHiiiE,)) =ay: s[El - {E, + Ej] y: 



(C25) 



(C26) 



To a good approximation the zero-phonon wave functions for J = 1 can be chosen to be 
composed of a single value of J^. Thus we may label the wave functions so that /i — m: 



{S^\{El)) = a E 5[El - {E, + EJ] 



J2cf,m{r,m)*Ci{r) 



The corresponding integrated intensity is 



7(°)^ f dE^is^ti) = ay: 



Ec/,m(r,m)*Q(r) 



(C27) 



(C28) 



The inner product of the J = wave function and the spatial part of the J — \ zero-phonon 
states is essentially unity. So 



(C29) 



h. Energy Gain by Conversion 

Here we give a similar analysis of the energy-gain spectrum at low temperature due to 
ortho-para conversion. The derivation is similar to that for para-ortho conversion so we only 
quote the results: 
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EQ,m(r,m)*c/(r) 



(C30) 



where is the probabihty that the J — 1, — m state is occupied and the role of initial 
and final states is interchanged from the para to ortho processes. The corresponding integrated 
intensity is 

2 



7(°)^ j dE{Sfl,) = BY.P„ 



EQ,m(r,m)*c/(r) 



= B . 



(C31) 



Since B = 3A we see that ratio of the total cross section for ortho to para conversion to that of 
para to ortho conversion is (1 — x) to x, where x is the ortho concentration. Normally the ratio 
of energy gain to energy loss cross sections follows the Boltzmann factor. Here, the populations 
are set by x rather than by the temperature. 



c. Single Phonon Energy Loss 
We have the powder average of 5*1^11.1 of Eq. (C21) as 

(^£i;i> = K^o E mE-E, + Ef)Y.\{f\ur.\i) 

Ji=l,Jf=l L 

and the corresponding integrated intensity is 



(C32) 



) = j dEr^{Sl%,iEr.)) = Ik'{u')Do . (CSS) 
In terms of the amplitudes of the wave function on the mesh points, the above result is 



(1) N 



■Do^PiS{EL-Ef + Ei)J2 



J2ci{M, r)c/(M, r)rL 



Also we obtain the powder average of S'^l^^.a and S'^j^^.g of Eq. (C21) as 
{S{%2{El)) = {Si%s{E^)r = -IJ^D.k' E -Ef + E,) 



M,r 



(CS4) 



X 



E C(112;M,M0(-l)^^+*^'(/|«_MTf+^'(J)K)(^|«_M'|/) , (CS5) 



MM' 



where the (J) are the operator equivalents of the spherical harmonics: 



SI 



7?'(J) = \/^4, T,^\:i)-T^^{J.J± + J±J.), T°(J) = ^^(3J,^-2). (C36) 
Here 

{i\u-M'\f) = X] ^^(a*' r)r_M' ■ (C37) 

and 

(/|«_Mrf+^'(J)|z) = ^XQ(/x,r)c^(Ai + M + M',r)V_MC(121;Ai,M + M') , (C38) 
where we used 

(J = 1; J, = M|T2^(p)| J = 1; J, = M') = -^6m,l+m'C{121; M', L) . (C39) 

VoTT 



We have the contributions to the integrated intensity 

--IMd^k'J^P, Y: C(112;M,M')(-l)^+^'(i|i.-M'|/)(/|^^-MTf+^'(J)|i) . (C40) 

Here the sum over final states should be restricted to single phonon states. Higher energy states 
make only a small contribution to this sum. So we make the closure approximation that the 

sum over |/) extends over all states, in which case 

= -lMD,n'J2P^ E C^(112;M,M')(-l)^(.|«l,,«-MTf +^'(J)K) . (C41) 

^ * -"-^ i M,M' 

As illustrated by Eq. (20), the initial state \i) is dominantly comprised of a single value of J^- 
Thus in Eq. C41 M + M' = dominates. In addition, the system is nearly isotropic. Then 
EmC'(112;M, -M)(-l)^|ii_Mp = 0. So, to a good approximation, 

4^) + = . (C42) 

We have made several approximations, but our result for the total integrated intensity will not 
be much affected by these approximations. 
Similarly, we get 
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327r 



(C44) 



-——02^ Yl PAEl -Ef + Ei)C{n2; M, M')C(222; M + M', iy){-l)^' 
2b\/21 ij,u,M,M' 

x{f\u^MT2-'{m{f\^M'T2-''-'''-^{JW ■ (C43) 
We now evaluate the integrated intensity 

= J dEL{S[%^{EL)) =h+t2, 

where 

h^^K'D2 y: mu-MT2-'{j)\f){f\uMmm{-i)''^' ■ (C45) 

Making the closure approximation this is 

1 nTT 
'"^ i,IJ.,M 

= lK'D2Y2Pi{i\^-MUM\i)i-l)^ = Ik'D2{u') . (C46) 

i 

Similarly 

3977- 

t2 = -i=D2K^ PiC {112; M,M')C {222; M + M',u) 

25v21 i,M,M'u 

x(-l)^+'^(i|t.^Tf +^'+^(J)x._MTr(J)N) . (C47) 

Again, we treat |i) as having a single value of J^, so that M + M' — 0. Also we again assume 
spatial isotropy, so that Em C{n2; M, -M){-1)^\u-m\^ = 0. Then t2 = 0. 

d. Zero Phonon Ortho Cross Section 

Finally wc consider the zero phonon contribution to the ortho to ortho cross section. Taking 
the powder average, we have 

= 47rc E E Ef) Y.im{mm-'{p)\f){-^r ■ (C48) 

Ji=l Jf=l n 

If 5 is the energy of the = ±1 levels relative to the = level and if the Boltzmann 
probability of the state \m > is Pm, then we may write 
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{Si\) = inCPoS{E + S) Y: 

m=±l L 



{m\Yi(p)m' + \{m\Yf\p)m 



which gives 



+A7rC5(E-5) ^ Pm 

m=±l 



{m\Y,\p)\0)\' + \{m\Yf\p)m 



(C49) 



(C50) 



The ratio of cross sections at energy gain to those of energy loss does satisfy detailed balance 
because within the species ( J = 1) we do maintain thermal equilibrium. 

3. INTENSITY RATIOS 



We develop an expression for the ratio, rp, defined to be the integrated intensity due 
to phonon creation divided by that due to para-ortho conversion. Using the results for /*^") 
obtained above we have 



rp 



{1 - x) J dEL{sSi,{Er^)) {i-x)m 

W{v?)D^ + Ik''{u'')D^ ( X 



\-x 



jo(|«:p)' + 2j2(|«:p)' 



x 



(C51) 



The ratio rc, defined to be the ortho to para conversion cross section in energy gain to that in 
energy loss due to para to ortho conversion, is given by 

X 



rc 



1-x 



(C52) 



Finally, we have rj=i, defined to be the total cross section (counting both energy loss and 
energy gain) for transitions within the J = 1 ground manifold divided by that due to para to 
ortho conversion, as 

X \ |C(Po + Pi) 



rj=i 



1-x 



3A 



(Po + Pi)rj=i(T = 0) , 



where 



(C53) 



(C54) 
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TABLES 

TABLE I. Effect of a spherical symmetric perturbation on U3 states. The Us states are charac- 
terized by N, the total number of harmonic phonons. Wave functions in a spherical potential are 
characterized by angular momentum K. Here we give the effect of the perturbation A{r/a)^ , where 
(r^) = 3a^ for the isotropic harmonic oscillator in three spatial dimensions. In the last column we 
give the shift in the average energy of the multiplet of states of a given value of N. 



N 


N,K 


Energy 


Avg E 





(0,0) 


15 A 


15 A 


1 


(1,1) 


35 A 


35 A 


2 


(2,0) 


75 A 


65 A 




(2.2) 


G3 A 




3 


(3,1) 


119 A 


105 A 




(3,3) 


99 A 
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TABLE II. Energy level systematics for a ( J = 0) molecule in an octahedral interstitial site of Cqq. 
Here we show the removal of degeneracy from a manifold of initial symmetry I to manifolds of final 
symmetry JT due to a perturbation V, as calculated in lowest-order perturbation theory. Here N is 
the total number of phonons, K is the angular momentum, and the other group theoretical labels are 
as in Fig. 1. We give a typical eigenfunction to illustrate the symmetry. Here d is the degeneracy 
(Dcg) of the manifold and = {x^) = (y^) = {z^)- Here the coordinate axes coincide with the cubic 
[100] directions. 



X 




V 






Deg 


Energy 


N = 2,K = 2 




k{x^ + y^ + z^- f r^) 




Eg 


2 


f KCT^ 








xy 


T2g 


3 




N = ?>,K = 3 






(x^ - |xr^) 


Tiu 


3 










x(y2 - 2;2) 


T2u 


3 










xyz 


Au 


1 


72 4 


N = 1,K = 1, Tiu 


rY^^in) 


X{xy + yz + zx) 




Eu 
Au 


2 
1 


2Aa2 


N = 2,K = 2, T2g 


xy 


\{xy + yz + zx) 




Eg 
Ag 


2 
1 


-Aa^ 
2Aa2 


N = 3,K = 3, T2u 




X{xy + yz + zx) 




Eu 


2 
1 


fAa^ 
-3Aa2 


N = 3,K = 3, Tiu 




\{xy + yz + zx) 




Eu 
Au 


2 
1 


AAa^ 
-fAa^ 


N = 3,K = 1, Ti, 


xir' - \a^) 


\{xy + yz + zx) 




Eu 
Au 


2 
1 


3A(t2 
-6Aa2 
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TABLE III. Phonon levels for a (J = 0) H2 molecule in an octahedral site for orientationally 
ordered and disordered C20 respectively. Energies (meV) are with respect to the ground state energy 
Eo,i. The symmetry of each manifold of degenerate levels can be read from Fig. 1. 









Octahedral 


Pa3 (Se) 


N 




a 






1 




1 


14.38 


13.16 


1 




2,3 


14.38 


14.47 


2 




1 


28.26 


26.69 


2 




2,3 


28.26 


27.49 


2 




4 


30.69 


30.41 


2 




5,6 


31.73 


31.39 


3 




1 


41.62 


40.10 


3 




2,3 


44.39 


42.40^ 


3 




4 


44.39 


42.42^ 


3 




5,6 


45.23 


44.40 


3 




7 


45.23 


45.83 


3 




8,9 


50.07 


49.36 


3 




10 


50.07 


49.85 



a) These energies are accidentally almost identical. However, group theory indicates that these 
levels are generically nondegenerate. 
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TABLE IV. Matrix elements of < A'^Q;|yl2(r)|01 > (in meV) for H2 in octahedral and Sg potential, 
respectively. Elements not listed are expected to be zero by symmetry. Numerically such elements 
were found to be very small. For this table the wave functions within each degenerate manifold were 
chosen to make the matrix elements of A2{r) as simple as possible. For Pa3 symmetry, the z-axis 
is taken to be the local three-fold axis. Therefore the octahedral wavefunctions are not necessarily 
identical to the Pa3 wavefunctions. 



< ArQ;|^5(r)|01 > Octahedral Pa3 (Se) 



N a 


T 






1 





(0,0) 


(-1.286,0) 


2 1 


U 


(0,0) 


(-0.506,0) 


2 1 


1 


(a,0) 


(0,0) 


2 1 


2 


(0,0) 


(0,0) 


2 2 


1 


(0,a) 


(-0.176,-0.106) 


2 3 


2 


(0,a) 


(2.413,-2.388) 


2 4 





(0,0) 


(-0.224,0) 


2 5 





(0,0) 


(0,0) 


2 5 


2 


m 


(-6.500,-1.754) 


2 6 







(0,0) 


2 6 


1 


(0,0) 


(-8.862,-1.400) 


a) Our numerical results give a 


= 0.859 and p = 5.662 


(in meV). 
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TABLE V. Non-zero matrix elements of < lajAJ > (in meV) for H2 in octahedral 

and Se potential, respectively. We also note that < la|^2('')|l/^ >=< l/3|^2('')|lo; > and 
< la|^^'^(r)|l/? >= (—1)"^ < lQ;|A2(r)|l/3 >*. For the octahedral symmetry, the wave functions 
within each degenerate manifold were chosen to make the matrix elements of A2{r) as simple as 
possible. 





< la\Al{r)\ip > 




Octahedral 


Pa3 (Se) 


T 


a 


(3 









1 


1 


(-3.853,0) 


(-2.020,0) 





2 


2 


(-3.853,0) 


(-1.326,0) 





3 


3 


( 7.708,0) 


(-1.326,0) 


1 


1 


2 


(0,0) 


(-2.109,-1.131) 


1 


1 


3 


(0,0.881) 


(1.137,-2.104) 


1 


2 


2 


(0,0) 


(1.534,2.440) 




2 


3 


(-0.881,0) 


(2.435,-1.533) 




3 


3 


(0,0) 


(-1.531, 2.442) 


2 


1 


1 


(-4.119,0) 


(0,0) 


2 


1 


2 


(0,0.880) 


(-0.225,-3.022) 


2 


1 


3 


(0,0) 


(-3.022,0.226) 


2 


2 


2 


(0,0) 


(0.560,0.835) 


2 


2 


3 


(0,0) 


(-0.830,0.556) 


2 


3 


3 


(0,0) 


(-0.565,-0.831) 
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TABLE VI. Shift of the center of gravity (CG) and sphtting (in nicV) of the {J = 1, N = 0) 
manifold when the nominally octahedral site has octahedral and Sq symmetry. 





Our Calculations 




Experiment 


Quantity 


Octahedral (Oh) 


Pa3 (Se) 


Pa3 (Se) 


Shift of CG ^ 


-0.134 


-0.141 


0.35 


Shift of CG ^ 




—0.16 


0.35 


Splitting 1st order'^ 





0.487 




2nd order^ 





-0.010 




totaF 





0.477 


0.70 


total^ 





0.46 


0.70 



a) Perturbation result of Eq. (16). 

b) Obtained by direct diagonalization of Eq. (3) 

c) Perturbation result of Eq. (17). 

d) Prom Ref. 5. 
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TABLE VII. Energy (in meV) of J = 1 single-phonon states where zero-of-energy is taken to be 2B. 



OclahcdrM] (Oh) 



Pa3 (S,0 





Full Mesh^ 


Perturbation^ 




Full Mesh^ 


Perturbation'' 


Tig 


13.11 


13.14 




12.59 


12.52 


Tig 


13.11 


13.14 


Eg 


12.82 


12.87 


Tig 


13.11 


13.14 


Eg 


12.82 


12.87 


T2g 


13.6 


13.68 


Eg 


13.47 


13.51 


T2g 


13.6 


13.68 


Eg 


13.47 


13.51 


T2g 


13.6 


13.68 


Ag 


14.20 


14.16 


^g 


15.61 


15.78 


Ag 


15.75 


15.31 


Eg 


16.46. 


16.60 


Eg 


16.60 


15.79 


Eg 


16.46. 


16.60 


Eg 


16.60 


15.79 



a) Solution to Eq. (3) for the three component wavefunction on a mesh of points. 

b) Solution to Eq. (5) using wavefunctions and energies for a J = molecule as previously 
determined numerically on a mesh of points. 
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TABLE VIII. Basis functions within the one phonon (J = 1) manifold. Here x, y, and z are the 
one phonon states with a single excitation in the phonon associated with the x, y, and z direction, 
respectively. In terms of the mj states (denoted \mj)) within (J = 1) we have X = (| — 1) — |l))/\/2, 
Y = + I - 1))/V2, and Z = |0). 



Oh Symmetry^ 



^2g 



■Ig 



(xZ + zX), (yZ + zY), (xY + yX) 
(xZ-zX), (yZ-zY), (xY-yX) 
{2zZ -yY -xX), (xX-yY) 
{xX + yY + zZ) 



Pa3 Symmetry*^ 



Eg and E| 
A, 



{xZ + zX, yZ + zY) , {xZ - zX, yZ - zY), (xX - yY, xY + yX) 
zZ, xX + yY, xY-yX 



a) The x, y, and z directions are taken to coincide with the four- fold axis of Oh. 

b) The z direction coincides with the three-fold axis of Se. 
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FIGURES 

SPHERICAL OCTAHEDRAL Pa3 (Sg) 




FIG. 1. Energy levels of a spherical (J = 0) II2 molecule confined in various ways. Here (g) is 
the degeneracy and the symmetry labels arc given. Left: The molecule is in a spherical harmonic 
potential V{r) = \k{x^ + + z^) or an anharmonic spherically symmetric potential (i. e. a generic 
spherically symmetric potential) . For a harmonic spherically symmetric potential the energy depends 
only on N, the total number of phonon excitations in the oscillators along the three coordinate 
directions. For a spherically symmetric potential eigenstates are characterized by their total orbital 
angular momentum L. Center: the molecule is in a potential appropriate to the octahedral interstitial 
site of orientationally disordered (Fm3m) solid Ceo- Right: the molecule is in a potential appropriate 
to the octahedral interstitial site of orientationally ordered (Pa3) solid Ceo, in which case the site 
symmetry is Se- The potentials used for the interstitial cases are discussed in Appendix A. 
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SPHERICAL SPHERICAL OCTAHEDRAL Pa3 (S5) 



NO ROTON- 

PHONON 

COUPLING 



E=4/3 r 



Eg: xX-yY, 



K=0 



;(2) 



2zZ-yY-xX ^2) 



Eg 



Ag: xX+yY+zZ.-- 

; XI) ' 



/(I) 



Ag 



J=l®N=l\ 



..•<5). 



E=2/i5 n 

K=2 \ 



■■■(3)- 



E=-2/3 r; 



K=l 



W Tig: xZ-zX, 
; \ yZ-zY, 
xY-yX 



\K3) 



■•■•(2) 



Ag 



Eg 



T2g: xZ+zX, \'\ gg 



yZ+zY, \(2) 



Ag 



E (meV) 
31 

30 
+ 29 



28 



■ 27 



xY+yX ■■•(1). 

FIG. 2. Removal of degeneracy as roton-phonon interactions are introduced and the site symmetry 
is lowered. The degeneracy is indicated by the number in parenthesis. At the far left is shown the 
completely degenerate level when spherical symmetry is assumed and no roton-phonon coupling is 
present. The next panel shows the effect of allowing roton-phonon interactions but preserving overall 
spherical symmetry. Here K is the total (orbital plus orientational) angular momentum. In the 
next panel spherical symmetry is lowered to octahedral symmetry which is appropriate for II2 in the 
octahedral interstitial site in orientationally disordered Cqo- The far right panel (and the energy scale) 
applies to the case of H2 in orientationally ordered Ceo in which case the site symmetry is Sg. 
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lxY + yX> 



K = 2,T 



2g 



IxX-y Y> 
K = 2, E„ 




FIG. 3. Translation-rotation wavefunctions for a (J = 1) H2 molecule in an octahedral intersti- 
tial site with one phonon. Here the plane of the paper is the x — y plane and for simplicity only 
the dependence in this plane is depicted. Each figure eight represents an \X) or \Y) orientational 
wavefunction and the sign associated with each lobe of this p-like function is indicated. Each orienta- 
tional wave function is multiplied by a translational wavefunction \y), or \z), where for instance 
\x) ~ xexp[— ^(x/a)^]. The presence of a phonon in the coordinate thus causes the wave function 
to be an odd function of ra, as one sees in the diagrams. As indicated in Fig. 2, the total angular 
momentum, K, which is the sum of the angular momentum of the phonon and that of rotation, is 
a good quantum number whose value is indicated. Top, Left: a K = 2, T2g function, top right: a 
K = 2, Eg function, bottom left: a K = 1, Tig function, and bottom right a = 0, Ag function. 
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E (meV) 



-28 



Tc / 



J=l N=l 



f 



B 




J=0N=1 ^ J=1N=0 
A / Ortho H2 (x) 



14 



I j=0 N=0 

Para H2(l-x) 



■-0 



FIG. 4. A schematic representation of possible transitions between the rotation-phonon energy 
levels that could be observed in a neutron scattering experiment. At low temperature, only the 
( J = 0, = 0) and ( J = 1, = 0) states are populated and therefore the transition Tc can not be 
observed at low temperature. The transition T^ is proportional to the coherent cross section of H2 
and therefore very small. The transitions T^ and T^ have comparable cross sections (see text for 
details). 
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FIG. 5. The calculated transition probabilities from the (J = 1,A^ = 0) levels to the nine-fold 
manifold of(J = l,A^ = l) levels at T= 4 K. For each pair of energies these transition probabilities 
represent the appropriate sum over degenerate levels. Note that there are at least eight transitions 
with comparable probability, suggesting that rich features could be observed in a neutron scattering 
experiment. 
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FIG. 6. The solid curve is the ratio rp of the single-phonon cross section to that from para-ortho 
conversion as given in Eq. (25) as a function of momentum transfer. The dotted Une shows 50 x r j=i 
(T=0) as given in Eq. (26). The experimental situation of Ref. 5 corresponds to a momentum transfer 
between 2 and 4 A^'^. 
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FIG. 7. Neutron energy loss spectrum (middle) at 4.2 K. The bottom curve is the result from our 
calculation using the WS77 potential. The dashed and gray lines are the contributions from rotational 
and vibrational excitations, respectively. The top curve is the spectrum after arbitrarily scaling Af^ 
by half, indicating that the orientational potential used in our calculations is too anisotropic as far as 
the center-of-mass motion is concerned. 
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2 6 10 14 18 

Neutron Energy Gain (meV) 



FIG. 8. Neutron energy loss spectrums obtained from various commonly used intermolecular po- 
tentials. For each potential we give the value of S, the splitting of the (J = 1) zero-phonon levels, 

which may be compared to the experimentally determined value 6 = 0.75 meV.^ Note that the average 
phonon energies of the first two potentials (top and middle curves) are way off from the observed value 
of 14 meV. 
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Neutron Energy Gain (meV) 



FIG. 9. Temperature dependent neutron energy gain spectrum of H2 in Ceo (the data is taken 
from Ref. 5). The inset shows ln(I//o) versus 1/kT, where / is the intensity of the feature at about 
28 meV. The slope of the Une indicates an activation energy barrier of 14.8 meV. 
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